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Abstract 



The production of charmed and beauty baryons in proton-proton collisions at 
high energies is analyzed within the modified quark-gluon string model. We 
present some predictions for the experiments on the forward beauty baryon 
production in pp collisions at LHC energies. This analysis allows us to find 
useful information on the Regge trajectories of the heavy (66) mesons and 
the sea beauty quark distributions in the proton. 
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\£l ' 1. Introduction 
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Various approaches of perturbative QCD including the next-to-leading or- 



der calculations (NLO QCD) have been applied to analyze the production of 
heavy flavour mesons and baryons in hard pp collisions, see for example jl|, |2j 
and references therein. However, it is impossible to use these approaches to 
study the hadron production at small values of the transfer because the run- 
ning QCD constant a s can be large. A rather successful description of various 
characteristics of heavy-quark hadroproduction processes can be obtained by 
using the approach for describing the soft hadron-nucleon, hadron-nucleus 
and nucleus-nucleus interactions at high energies based on the topological 
1/N expansion in QCD [3|, |4|, where N is the number of flavours or colours, 
and the closely related string and colored-tube models [5|, |6|. There are suc- 
cessful phenomenological approaches for describing the soft hadron-nucleon, 
hadron-nucleus and nucleus-nucleus interactions at high energies based on 
the Regge theory and the 1/N expansion in QCD, for example the quark- 
gluon string model (QGSM) [7J, the VENUS model |8|, the dual parton model 
(DPM) 0. 
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The main components of the QGSM and the DPM are the quark distribu- 
tions in a hadron and the the fragmentation functions (FF) describing quark 
fragmentation into hadrons. These are expressed in terms of intercepts of 
the Regge trajectories or(0). The largest uncertainty in the calculations of 
the cross sections for the yields of heavy flavours in these models is mainly 
due to the absence of any reliable information on the transfer t dependence 
of the Regge trajectories of heavy quarkonia (QQ). Assuming linearity of 
the (QQ) trajectories, the intercepts turn out to be low, a^(0) = —2.2, 
Q!t(0) = —8, —16, and so the contribution of the peripheral mechanism de- 
creases very rapidly with increasing quark mass, Accordingly, the finding of 
the ^-dependence for oi(QQ)(t) in the region < t < M? q q, and estimations 
of their intercepts become especially important. 

In this paper we investigate the open charm and beauty baryon produc- 
tion in pp collisions at LHC energies and small pt within the QGSM to find 
the information on the Regge trajectories of the heavy (cc) and (bb) mesons, 
the FF of all the quarks and diquarks to heavy flavour hadrons and the sea 
beauty quark distribution in the proton. We also present some predictions 
for the future LHC experiments on the forward A^ production in pp collisions. 

2. General formalism for hadron production in pp collisions within 
QGSM 

Let us briefly present the scheme of the analysis of the hadron production 
in the pp collisions within the QGSM including the transverse motion of 



quarks and diquarks in colliding protons [10|]. As is known, the cylinder type 
graphs for the pp collision presented in Fig.l make the main contribution to 
this process J7|. The left diagram of Fig.l, the so-called one-cylinder graph, 
corresponds to the case where two colourless strings are formed between 
the quark/diquark (q/qq) and the diquark/quark (qq/q) in colliding protons; 
then, after their breakup, qq pairs are created and fragmentated to a hadron, 
for example, D meson. The right diagram of Fig.l, the so-called multicylinder 
graph, corresponds to creation of the same two colourless strings and many 
strings between sea quarks/antiquarks q/q and sea antiquarks/quarks q/q in 
the colliding protons. The general form for the invariant inclusive hadron 
spectrum within the QGSM is @, El} 

E ^=^^r^ an{s)Mx ' Pth (1) 
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Figure 1: The one-cylinder graph (left) and the multicylinder graph (right) for the inclusive 
pp — > hX process. 



where E, p are the energy and the three-momentum of the produced hadron 
h in the laboratory system (l.s.); E*,s are the energy of h and the square 
of the initial energy in the c.m.s of pp; x,pt are the Feynman variable and 
the transverse momentum of h; a n is the cross section for production of the 
n-Pomeron chain (or 2n quark-antiquark strings) decaying into hadrons, cal- 
culated within the quasi- "eikonal approximation" |ll| . Actually, the function 
<fi n (x,pt) is the convolution of the quark (diquark) distributions and the FF, 
see the details in |7| and [9|, ll(| . 



€ P (^) = FW(x + )FW(x_) + FW(x + )FW(z_) + 2(n - l)^^)^^.) , (2) 
where x± = \ ( \J x 2 + x\ ± x) , and 

F^ n \x ± ) = j dx l fi n \x 1 )G T ^ h (j^j , (3) 

Here r means the flavour of the valence (or sea) quark or diquark, fr (%i) 
is the quark distribution function depending on the longitudinal momen- 
tum fraction x\ in the n-Pomeron chain; G T ^h(z) = zD T ^ h (z), D T ^ h (z) 
is the FF of a quark (antiquark) or diquark of flavour r into a hadron h 
(charmed hadron in our case); o n is the cross section for production of the 
n-Pomeron chain (or 2n quark-antiquark strings) decaying into hadrons, cal- 
culated within the "eikonal approximation" |7|. 

To calculate the interaction function (f) n (x,p t ) we have to know all the 
quark (diquak) distribution functions in the nth Pomeron chain and the 
FF. They are constructed within the QGSM using the knowledge of the 
secondary Regge trajectories, see details in [7|, [l0|. Note that the modified 
QGSM including the internal longitudinal and transverse motion of quarks 



(diquarks) in the initial nucleon has been successfully applied [12|]-[14| to 
describe the production of heavy flavour hadrons produced in p — p and p — p 
collisions at not large pt- 

3. Heavy baryon production within QGSM 

3.1. Sea charm and beauty quark distribution in proton 



According to the QGSM, see for example [12| and references therein, the 
distribution of c(c) quarks relevant to the nth Pomeron chain (Fig. 1 (right)) 
is, see for example [l2| and references therein, 

f^)=Cip ci e>X^(l-xY- (4) 

where a cn = -cty(O), g cn = 2(a p (0) - a B (0)) - aty(O)) + n - 1; S c{d ) is the 
weight of charmed pairs in the quark sea, Clfl is the normalization coefficient 



'0(c) 



10J, 0^(0) is the intercept of the ip- Regge trajectory. Its value can be 
—2.18 assuming that this trajectory a^(t) is linear and the intercept and the 
slope ct4(0) can be determined by drawing the trajectory through the J/^/- 
meson mass mj/q, ~ 3.1 GeV and the x-meson mass m x = 3.554 GeV [151 ] . 
Assuming that the ip- Regge trajectory is nonlinear, one can get a^(0) ~ 0, 
which follows from the QCD-based model analysis [la, |l7|. The distribution 
of b(b) quarks in the nth Pomeron chain (Fig. 1 (right)) has the similar form 

where a bn = -a T (0), g bn = 2(a p (0) - a B {0)) - a T (0)) + n - 1; a p (0) = 
1/2 is the well-known intercept of the p-trajectory; ob(0) — —0.5 is the 
intercept of the baryon trajectory, a-r(O)) is the intercept of the T- Regge 
trajectory, its value also has an uncertainty. Assuming its linearity, one can 



get «t(0)) = —8,-16 [21], while for the nonlinear (bb) Regge trajectory 
«x(0) — [l6|, [l7|. Inserting these values in the form for f<?l(x) and f b ^(x) 
we get the large sensitivity for the c and b sea quark distributions in the nth 
Pomeron chain. Note that the FFs also depend on the parameters of these 
Regge trajectories. Therefore, the knowledge of the intercepts and slopes 
of the heavy-meson Regge trajectories is very important for the theoretical 
analysis of open charm and beauty production in hadron processes. 

Note that all the quark distributions obtained within the QGSM are dif- 
ferent from the parton distributions obtained within the perturbative QCD 



which are usually compared with the experimental data on the deep inelas- 
tic lepton scattering (DIS) off protons. To match these two kinds of quark 
distributions one can apply the procedure suggested in [l8(. The quantities 
g cn or gi, n entering into Eq.(TJ| and Eq.(j5]) are replaced by the following new 
quantities depending on Q 2 

~g cn = g cn (i + ^) ; hn = g bn (l + ^) (6) 

The parameters c and d are chosen such that the structure function con- 
structed from the valence and sea quark (antiquark) distributions in the 
proton should be the same as the one at the initial conditions at Q 2 = Qq 
for the perturbative QCD evolution. A similar procedure can be used to get 



the Q 2 dependence for the powers a cn and o^ entering into Eqs.( H|5]) [18 



Then, using the DGLAP evolution equation [19j, we can obtain the structure 
functions at large Q 2 . 

3.2. Fragmentation functions of partons to heavy flavour baryons 



According to [10] . the production of B- baryons, for example A, A c , A&, is 
illustrated in Fig.2. Fig. 2a illustrates the B production in the fragmentation 
region, when the diquark in the colliding proton picks up directly the b quark 
from the qq chain created between the diquark in the first colliding proton 
and the quark in the second proton. Fig. 2b also illustrates the production 
of B in the fragmentation region, when one quark and the string junction in 
the first proton pick up a quark (u or d) and the b quark from the qq chain 
mentioned above. Figure. 2c corresponds to the baryon-antibaryon (BB) pair 
production in the central region. 

The FF of quark (diquark ) c hains to heavy flavor hadrons were con- 
structed within the QGSM in J20j, l2l| . For example, the FF of quarks (an- 
tiquaries) to the charm and bottom baryons A c , A^ have the following forms 
0H : 



Q A c( A b)( z \ — G Ac( - Ab \z) = a Ac(Ab) (l - /2 ) _Q V'(T)(0)+A+2(a i? (0)-Q JV (0)) /y\ 

and 

G AaiAb \z) = Gf iAb) (z) = G A ^ = Gt (Ab) = G A ^ A »\z)(l - z)^- a «^ (8) 



c c c 



a) b) c) 

Figure 2: Diagrams corresponding to the fragmentation of diquarks into baryons 
-B(A,A C ,A(,) in the fragmentation region (a) and in the central region (b), when the 
baryon-antibaryon pair (BB) is produced, the dashed line corresponds to the string junc- 
tion between the diquark and the quark. 



The FF of diquarks to A c .Aj, are the sums of two parts, central G Qqq ) and 



fragmentation G / . These parts have the following forms 22 



QA c (A b ) = a M^b)^ _ z \-a^, m (Q)+\-4[l-a N (0)} /g\ 

and 

Gf u [ Ab) = af {Ab) z 2 (l - z )-^m(o)+x+x ; (10) 

£,A c (A b ) _ ^A c (A b ) _ a A c (A 6 )^_ 2 ( Q!JV (0)-a i? (0)/ 1 _ ^-^ (T )(0)+A q-q 

where A = 2o/(0) R < pj >~ 0.5. To calculate inclusive spectra of A c or A b 
baryons we have to know the constants a and a* . The constants 



a c and a, c are presented in [21], their application allowed describing the 
experimental data on the A c production. However, the similar constants for 
the At production are unknown. The question arises how to find these values. 



The constant a f c at a^,(0) = was taken in [2l| as the same value as for the 
A hyperon production a A [lOj, e.g., a Ac = a A = 0.02 assuming that in the 
fragmentation region (Fig.2(a,b)) the production mechanism for A is the same 
as for the A c production. We can assume it for the A& production as well, e.g., 
aj b = a f c = a A = 0.02. As for the A& production in the central region, when 

the A fe -A fe pair is produced, see Fig. 2c, one can estimate a Q b from the relation 



of the A& production in the central region to the A c production. In the central 
region (Fig. 2c), in addition to the production of two pairs of light quarks- 
antiquarks (u, d), the pairs of heavy quarks- ant iquarks ((ss), (cc), (bb)) are 
produced, when there are (AA), (A C A C ) and (A^A;,) pairs in the final state. 
Knowing the probability w qq for the (qq) pair production from the decay of 
the string between the diquark and the quark (Fig.2) one can relate a c to 
a A and a b to a c , e.g., 

4° = ^4 (12) 



and 



!M a A. = !M a A (13) 

Wrr. W. q z 



The probability of the qq pair production per unit of the four- dimension 
volume of the color tube w qq was calculated in [6[ within the colour-tube 
fragmentation model p, la] in the following form: 

w m = T^ e ^( _7rm g/ 'P) ' ( 14 ) 



where p = l/{2 r wa' R ) is the so-called string tension, a' R is the slope of the 
Regge trajectory, m q is the quark mass. Note that Eq.(14) was obtained in 
[6] on the assumption of the linearity of the Regge trajectories. Calculating 
w c5 and w sS using Eq. (IT4"|) and taking the value a A = 0.3 from 



10 



one can 
get a c = 4 ■ 10 -4 that coincides with the value presented in 2l| which was 
used both at a*(0) = —2.18 and a#(0) = to describe rather successfully 
the experimental data on the A c production in pp collisions. It makes us 
confident that one can use the same method to estimate a b as well. Taking 
a'j(0) ~ 0.1 (GeV/cp 2 from 2l[ and calculating Eq. (I13J14J) one can get 
a b ~ 1.7 • 10 -6 , which we also use both at or(0) = —8 and a-f(0) = 0. 

3.3. Charmed and beauty baryon production in pp collision 

The information on the charmonium (cc) and bottomonium (bb) Regge 
trajectories can be found from the experimental data on the charmed and 
beauty baryon production in pp collisions at high energies. For example, 
Fig. 3 illustrates the sensitivity of the inclusive spectrum da/dx of the pro- 
duced charmed baryons A c to different values for a^(0). The solid line cor- 
responds to aty(0) = 0, and the dashed curve corresponds to a^(0) = —2.18. 



Measuring the decay A c — > A°37r, the i?608 experiment 23] at | x |> 0.5 and 
y/s = 62 GeV obtained the pp -> A C X cross section 2.84±0.50±0.72/jb. The 
branching ratio of the decay A c — > A°3tt is 2.8±0.7±1.1% therefore, the cross 
section of the A c production is cr(| x |> 0.5) = 101 ± 18 ± 26/zb. Theoretical 
expectation for this cross section is <r(| x |> 0.5) = 87.3/zb at a$(0) = and 
cr(| x |> 0.5) = 30.5/do at a*(0) = -2.18. On the other hand, the R422 
experiment [24| was measuring the cross section of the process p — > e~A c X 
at | x |> 0.35 with a large uncertainty from (26 ± 12) fib to (225 ± 9) fib 
depending on the assumption made, see Table 6 in [24j. So, the open circles 
in Fig. 3 corresponding to the R608 experiment are more adequate to our 
calculations at a^(0) = 0, see the solid line in Fig.3. 
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Figure 3: The differential cross section da/dx for the inclusive process pp — > A C X at 
\fs = 62 GeV. The solid line corresponds to a^(0) — 0, and the dashed curve corresponds 
to Oty(0) = —2.18. The open circles correspond to the R608 experiment (23j, and the dark 
circles correspond to the R422 experiment [24j | . 



High sensitivity of the inclusive spectrum da/dx of the produced beauty 
baryons A& to different values for «x(0) is presented in Fig. 4 (left). The p t - 
inclusive spectrum of A& has much lower sensitivity to this quantity, according 
to the results presented in Fig.4 (right). Actually, our results presented in 
Fig.4 could be considered as some predictions for future experiments at LHC. 

The produced A& baryon can decay A5 -> J/\&A° with the branching 
ratio Br Ab ^j/y A o = T Ab ^ J/q , A o/T tot = (4.7 ± 2.8) ■ 10 -4 and J/* decays 
into fi + fi~ {Brj/xT f ^. fi + fi _ = (5.93 ± 0.06)%) or into e + e~ (Brj/y^. e + e - = 



5.93 ±0. 



whereas A can decay into pn (Br A o 



—>piT~ 



63.9 ±05%), or 
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Figure 4: The differential cross section da/dx (left) and da/dP? (right) for the inclusive 
process pp — > Ai,X at y/s = 10 TeV. 



into mi (5r ao- 



35.8 ± 0.5%). Experimentally one can measure the 



differential cross section da/d^ p dt p dMj/^, where £ p = Ap/p is the energy- 
loss, t p = (p in — pi) 2 is the four- momentum transfer, Mj/y is the effective 
mass of the J/\l/-meson. 



4. Results and discussions 

Now let us analyze the production of the beauty hyperon, namely A°, 
at small scattering angles 9\ h in the pp collision at LHC energies. This 
study would be reliable for the forward experiments at the LHC. In Fig.5 
the distributions over 9j/y (left) and 9^+ (right) are presented at different 
values of the intercept ax(0) = (solid line) and ax(0) = —8 (dashed line), 
where 9j/^ is the scattering angle for the final J/\l/. Figure 5 shows a large 
sensitivity of these distributions to the intercept «x = and ax = —8 of the 
T(bb) Regge trajectory at small values of the scattering angle. It also shows 
that the angle distribution of the \x meson produced from the J/\l/ decay is 
broader than the one for the J/\l/ meson. 

In Fig. 6 the distributions over 9\ b (left) and 9 P (right) are presented at 
different values of the intercept ar(0) = (solid line) and a~f(0) = —8 
(dashed line), where 9\ b is the scattering angle for the A& baryon and 9 P is 
the scattering angle of the final proton. One can see from Fig. 6 that the 
angle distribution of protons produced from the A decay is broader than 
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Figure 5: The distributions over 



v/* 



(left) and # M + (right) for the inclusive process 



pp —> AbX — >• fi + fi pn X at y/s — 10 TeV, where 6j/^ and 8^+ are the scattering angles 
for J/* and fi + . 



the one for the A hyperon. In Fig. 7 the same 6 P distribution as in Fig. 6 is 
presented but at very low values, 6 „ < 1.6 mrad, which can be observed in 
the TOTEM experiment at CERN |26 



All the angle distributions presented in Figs. 5-7 show a huge sensitivity 
to the values of the intercepts of the T(bb) Regge trajectory «t(0) = and 
a-f (0) = —8 at the angles about 4 — 8 mrad for muons, pions and protons 
and at the angles about 2 — 4 mrad. for A& baryons and J/\l/ vector mesons 
with the hidden charm. The distribution over the variable £ p related to the 
proton loss energy presented in Fig. 7 (right) shows that the cross section 
for «x(0) = —8 is too low at C, p < 0.6 and, in fact, it cannot be measured 
experimentally. 

From the distributions in Figs. 5-7, the following experimental signatures 
can be deduced. 

The ATLAS forward detectors could register the decay A° — > J/\& A — > 
[i + fi~ TT°n by detecting the neutron in the ZDC 25) and two muons in the 
muon detector. This possibility is under study. 

The TOTEM [26| together with the CMS might be able to measure the 
channel A^ — > J/^S A — > fi + {i~ 7i~p (the integrated cross-section is about 
0.2-0.3 fib at ct T (0) = and smaller at a T (0) = -8). The T2 and Tl 
tracking stations of the TOTEM apparatus have their angular acceptance 
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Figure 6: The distributions over 6^ b (left) and 
A b X -> fi+^-pir-X at y/s = 10 TeV. 



(right) for the inclusive process pp 



in the intervals 3 mrad < 9 < 10 mrad (corresponding to 6.5 > i] > 5.3) 
and 18 mrad < 9 < 90 mrad (corresponding to 4.7 > rj > 3.1) respectively, 
and could thus detect 42% of the muons from the J/\l/ decay (Fig. 5). In 
the same angular intervals, 36% of the n~ and 35% of the protons from the 
A decay are expected. According to a very preliminary estimate [271 ]. pro- 
tons with energies above 3.4 TeV emitted at angles smaller than 0.6 mrad 



could be detected in the Roman Pot station at 147 m from IP5 26|, |27| . In 
the latter case, the reconstruction of the proton kinematics may be possible, 
whereas the trackers Tl and T2 do not provide any momentum or energy 
information. Future detailed studies are to establish the full event topolo- 
gies with all correlations between the observables in order to assess whether 
the signal events can be identified and separated from backgrounds. These 
investigations should also include the CMS calorimeters HF and CASTOR 
which cover the same angular ranges as Tl and T2 respectively [271 ]. 



5. Conclusion 

We analyzed the production of charmed and beauty baryons in proton- 
proton collisions at high energies within soft QCD, namely the quark-gluon 
string model (QGSM). This approach can describe rather satisfactorily the 



charmed baryon production in pp collisions [12j-ll4j. It allows us to apply the 
QGSM to studying the beauty baryon production in pp collisions. We focus 
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Figure 7: The distributions over 9 p for the inclusive process pp — > A b X — > fi + /j,~ pn~ X at 
^s = 10 TeV. The solid curve corresponds to a T (0) = 0; the long dashed line corresponds 
to a-r(O) = —8 (left). Right distributions over £ p = (y / i — 2E p )/^/s for the inclusive 
process pp — > AbX — > (i + [i~pTr~X at ^/s — 10. TeV. 



mainly on the analysis of the forward A& production in pp collisions at LHC 
energies and got some predictions which could be reliable at the TOTEM 
and ATLAS experiments at CERN. It is shown that information on the sea 
6-quark distributions in the proton and the fragmentation functions of quarks 
(diquarks) into beauty baryon A& can be extracted from this analysis. We also 
show that the angle distributions of the hadrons produced in the reactions 
pp — > A b X — > n + n~pii~X or pp — > A b X — > fi + fi~mr°X calculated within 
the QGSM are very sensitive to the intercept values «x(0) = 0, —8 of the 
T(bb)- Regge trajectory. It means that the experiments on the forward A^ 
production in pp collisions at the LHC could give us the answer whether the 
T(bb)- Regge trajectory is linear (a-r(O) = —8) or nonlinear (a-r(0) = 0); 
the latter is predicted by the perturbative QCD. We did not include the 
diffractive and double diffractive A& production in pp collisions because these 
pro cesses can be experimentally separated from the forward A^ production 

Note that in this paper we neglect the contribution of the intrinsic charm 
in the proton calculating the charmed baryon production in pp collisions and 
possible intrinsic beauty in the proton. However, as shown recently 28| . 
the intrinsic charm in the proton can result in a sizable contribution to the 
forward charmed meson production. Therefore, we intend to include this 



12 



effect in the next more detailed study of this problem. 
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